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THEORETICAL  AERODYNAMIC  CHARACTERISTICS  OF  A  FAMILY 
OF  SLENDER  WING -TAIL -BODY  COMBINATIONS 
By  Harvard  Lomax  and  Paul  F .  Byrd 


SUMMARY 


The  aerodynamic  characteristics  of  an  airplane  configuration 
composed  of  a  swept-back  wing  and  a  triangular  tail  mounted  on  a 
cylindrical  body  are  presented.  For  simplicity,  the  leading  edge  of 
the  wing  is  considered  to  be  straight  and  the  trailing  edge  to  be 
shaped  so  that  the  span-loading  curve  is  flat  between  the  fuselage  and 
the  wing-tip  regions;  the  result  is  a  nearly  constant-chord  swept-back 
wing.  A  method  by  which  other  trailing-edge  shapes  can  be  studied  is 
indicated.  The  analysis  is  based  on  the  assumption  that  the  free-stream 
Mach  number  is  near  unity  or  that  the  configuration  is  slender.  The 
calculations  for  the  tail  are  made  on  the  assumption  that  the  vortex 
system  trailing  back  from  the  wing  is  either  a  sheet  lying  entirely  in 
the  plane  of  the  flat  tail  surface  or  has  completely  "polled  up"  into 
two  point  vortices  that  lie  either  in,  above,  or  below  the  plane  of  the 
tail  surface. 


INTRODUCTION 


The  study  of  lifting  surfaces  flying  at  either  subsonic  or  super¬ 
sonic  speeds  at  small  angles  of  attack  has  been  reduced,  by  the  well- 
known  process  of  linearization,  to  the  study  of  the  equation 

(1-Mo2)  cpxx  +  9yy  +  Tzz  =  0  (l) 

where  9  is  a  perturbation  velocity  potential  in  a  field  having  a 
uniform  free-stream  velocity  V0  directed  parallel  to  the  x  axis, 
and  where  Mo  is  the  Mach  number  of  the  free  stream. 

One  basic  simplification  of  equation  (l)  is  brought  about  by 
neglecting  velocity  gradients  along  the  span  of  the  wing.  If  the  wing 
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is  lying  in  the  z  =  0  plane,  this  amounts  to  neglecting  the  term  cpyy 
in  equation  (l),  and  results  in  the  'well-known  partial  differential 
equation  by  means  of  which  two-dimensional  or  section  characteristics 
are  studied. 

Another  basic  simplification  of  equation  (l)  can  be  attained  by 
neglecting  the  term  (l-M02)  Txx*  Such  a  procedure  is  possible  when  the 
Mach  number  is  close  to  1  or  the  wing  plan  form  is  so  slender  that 
velocity  gradients  in  the  free-stream  direction  are  negligible  in  com¬ 
parison  with  the  gradients  in  the  y  and  z  directions.  Equation  (l) 
has  already  been  analyzed  in  these  two  connections  in  references  1  and  2 
for  certain  plan  forms.  The  purpose  of  this  report  is  to  extend  this 
theory,  which  has  been  named  slender  wing  theory,  to  include  an  entire 
airplane  configuration. 

Results  are  presented  for  a  nearly  constant-chord,  swept-back  wing 
mounted  on  a  cylindrical  body  having  a  triangular  horizontal  tail  located 
aft  of  the  wing  trailing  edge.  Both  wing  and  tail  are  flat  surfaces, 
and  the  results  are  only  those  due  to  changes  in  the  airplane  angle  of 
attack. 

A  list  of  important  symbols  is  given  in  appendix  A. 


I  -  SWEPT-BACK  WING  ON  A  BODY  OF  REVOLUTION 

Partial  Differential  Equation,  Boundary 
Conditions,  and  Form  of  the  Solution 


Under  the  assumption  that  the  free-stream  Mach  number  is  1  or  that 
the  perturbation  velocity  gradient  in  the  x  direction  is  small,  the 
partial  differential  equation  which  must  be  satisfied  for  the  solution 
of  lifting  surface  problems  can  be  written 

Tyy  +  9zz  =  0  (2) 

Equation  (2)  is  simply  Laplace’s  equation  in  two  dimensions,  the 
variables  representing  lateral  and  vertical  coordinates  in  a  plane 
transverse  to  the  direction  of  motion. 

The  boundary  conditions  associated  with  equation  (2)  are  given 
along  a  line  in  this  transverse  plane  and  specify  that  the  fluid  veloc¬ 
ity  is  everywhere  tangential  to  the  surface  of  the  body.  The  problem  is, 
of  course,  to  find  at  other  points  in  the  plane  the  potential  that 
satisfies  equation  (2)  and  fits  these  boundary  conditions.  Of  particular 
interest  is  the  streamwise  component  of  velocity  along  the  surface  of 
the  wing  and  body  since  this  is  directly  related  to  the  loading  thereon. 
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Solutions  to  equation  (2)  are  readily  available.  Two  different 
analytic  forms  of  these  solutions  will  be  used  in  the  following  analysis 
One  form  is  concerned  with  the  use  of  the  complex  variable,  the  other 
with  the  use  of  Green’s  theorem  and  the  inversion  of  a  real,  singular, 
integral  equation.  In  general  the  procedure  will  be  to  use  concepts 
associated  with  the  complex  variable  to  map  the  boundary  conditions  onto 
a  slit  along  the  real  axis,  then  to  solve  the  resulting  problem  by 
inverting  an  integral  equation,  and  finally,  to  use  the  complex  variable 
again  to  extend  such  a  solution  out  into  space  by  the  principle  of  ana¬ 
lytic  continuation. 


Discussion  of  notation  and  transformations.-  The  first  part  of  this 
report  will  be  devoted  to  the  analysis  of  the  configuration  shown  in 
sketch  (a).  The  following  is  a  description  of 
this  configuration.  Everywhere  behind  the 
leading-edge-fuselage  juncture  the  fuselage 
is  a  circular  cylinder  having  a  radius  rQ. 

Ahead  of  this  juncture  the  fuselage  comes  to 
a  point,  the  manner  being  arbitrary.  The 
wing  is  a  flat  plate  without  twist  or  camber 
mounted  at  zero  incidence  on  the  fuselage 
and  the  whole  configuration  is  placed  at  a 
small  angle  of  attack  a  with  respect  to 
the  free-stream  direction.  The  origin  of 
the  coordinate  system  is  located  at  the 
wing  apex.  The  leading  edge  is  a  straight 
line  with  slope,  dy/dx,  equal  to  m.  It 
will  be  convenient  at  some  places  in  the 
report,  however,  to  use  the  expression 
y  =  s(x)  for  the  equation  of  the  leading  (a) 

edge,  hence,  s(x)  and  mx  are  used  inter¬ 
changeably.  The  trailing  edge  is  repre¬ 
sented  by  the  line  y  =  t(x)  and  is,  in 

general,  not  straight.1  The  maximum  semispan  of  the  wing  is  denoted 
by  sQ.  The  symbol  tQ,  as  can  be  seen  in  the  sketch, ' refers  to  the 
lateral  distance  from  the  x  axis  to  the  point  at  which  the  trailing 
edge  intersects  a  line  that  is  parallel  to  the  y  axis  and  passes 
through  the  last  outboard  point  of  the  leading  edge.  Finally,  c0  is 
the  chordwise  distance  from  the  origin  to  the  trailing-edge -fuselage 
juncture . 


i-It  was  considered  advisable  at  this  time  to  consider  only  the  rather 
particular  configuration  outlined.  As  the  analysis  progresses  it 
will  be  pointed  out  where  the  solution  can  be  generalized  to  include, 
for  example,  wings  with  twist  and  camber. 
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A  second  coordinate  system  will  also  be  used  in  the  succeeding 
development.  Let  the  y,z  plane  be  represented  by  the  complex 
variable  g , 

g  =  y  +  iz  =  pe10 
then  introduce  the  plane, 

61  =  Vl  +  izi  = 

so  that  the  £  plane  maps  onto  the  plane  by  means  of  the  Joukowski 

transformation 

lx  - 1  +  rf  (3) 

By  means  of  such  a  transformation,  the  circle  of  radius  rQ  which 
represents  a  section  of  the  fuselage  in  the  I  plane  maps  onto  a  portion 
of  the  real  axis  in  the  |i  plane  (see  sketch  (b)); 


£  plane 


and  the  part  of  the  real  axis  which 
lies  outside  the  circle2  in  the  i 
plane  maps  into  the  remaining  part 
of  the  real  axis  in  the  plane. 

As  a  consequence  of  equation  (3) 7 
the  following  relations  hold  for  zx 
equal  to  zero: 


^The  Joukowski  transformation  is  double  valued  in  that  the  regions 
inside  and  outside  the  circle  p  =  Tq  both  map  onto  the  entire  1 1 
plane.  In  this  report  only  the  field  outside  the  circle  is  of  interest. 
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and 


Further , 


yi  =  y  +  r-~>  yi2>ri2 

yi  =  2r0  cos  6,  yi2<ri2 


r  r 

Si  =  S  +  tx  =  t  +  rx  =  2r0 


(5) 


(6) 


From  the  basic  theory  underlying  the  use  of  complex  variables  in  fluid- 
flow  theory,  induced  velocities  in  the  two  planes  are  related  by  the 
expression 


V-1W  =  (  V 1  -  1W  )  - = 

d£ 


(7) 


from  which,  since  in  polar  coordinates 


dll 

d| 


cos  20 


+  1 


sin  20 


(8) 


it  follows  that 


V  =  Vi 


w  =  w. 


K?)2 


cos  20 


1  -  m 


cos  20 


<rj 


sin  29 


sin  26 


Vy  = 


)  (9) 


=  vx  cos  0+W!  sin  9  1-  ^2^  +2wx  sin  0 


Lastly,  Laplace’s  equation  must  also  be  satisfied  in  the  ii  plane, 
hence 


('Pihxzi  +  0 


(10) 
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Boundary  conditions.-  In  this  part  of  the  report  (part  I),  the 
effect  of  a  cylindrical  body  mounted  on  a  nearly  constant-chord  swept- 
back  wing  will  be  studied.  (Reference  3  contains  an  analysis  of  the 
effect  of  a  cylindrical  body  mounted  on  a  triangular  wing,  and  reference  Ij- 
presents  results  for  a  swept-back  wing  with  no  body;  both  references 
use  the  assumptions  of  slender -wing  theory.)  The  boundary  conditions 


will  be  presented  in 

the  y,z 

space  first 

and  can  be  written 

(i) 

0 

11 

> 

P  = 

r0,  0<9<2n 

(ii) 

0 

11 

> 

z  = 

p  P  P 

0,  t  <  y  <  s  > 

(11) 

(iii) 

<1 

11 

0 

X. 

w  =  V0a,  P  = 

«  ,  0  <  0  <  2n 

Equations  (ll)  represent  the  conditions  for  a  cylinder  located  at 
p  =  r0  and  two  wing  panels  located  between  ±t  and  ±s  on  the  real 
axis,  both  cylinder  and  wings  being  at  rest  in  a  free  stream  which  is 
moving  with  velocity  w  =  V0a  at  infinity. 

It  follows  from  equations  (9)  that  these  boundary  conditions 


become,  in  the  plane, 

(i) 

3 

H 

II 

O 

>• 

2  2 
0<yi  <  ri 

di) 

O 

II 

H 

2  2  2 

ti  <  yx  <  si 

'  (12) 

(iii) 

▼1  =  0,  w-L  =  vca, 

Pi  =  ® 

li 

Equations  (12)  represent  the  boundary  conditions  for  three  wing  panels 
along  the  real  axis,  all  at  rest  with  respect  to  the  free  stream  moving 
with  velocity  wx  =  VQa  at  infinity  in  the  transformed  plane. 


It  is  more  convenient  to  work  with  boundary  conditions  which  vanish 
at  infinity,  however,  so  the  final  form  of  the  conditions  which  must  be 
satisfied  is  derived  from  equations  (12)  by  subtracting  the  free-stream 
velocity  V0a.  There  results 

(i)  vx  =  -  VDa,  0<y12<r12 

(ii)  wx  =  -  VQa,  t12<y12<s12 

(iii)  vx  =  wx  =  0,  px  =«»  ,  0 <9<2k 

t 

General  solution.-  The  general  solution  to  equation  (10)  which 
gives  the  vertical  induced  velocity  wx  at  a  point  in  the  lx  plane 
due  to  the  jump  in  the  value  of  the  induced  velocity  v±  across  the  y1 
axis  can  be  written  (see,  e.g.,  reference  5) 
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wi(yi>zi) 


1_  n*1  (yi-ya)^vi(ya) 

2ltV  (yx-yaf+Zi2 

-Si  “■  ± 


<3-y2 


(14) 


where  y2  is  the  variable  of  integration.  Set  zx  equal  to  zero  and 
there  results  the  value  of  the  vertical  induced  velocity  on  the  y1 
axis.  Thus 

”i(yi>  =  wi(yi> °>  ■  -  tf  ^  (15) 

"S1 

Equation  (15)  is  the  form  of  the  solution  which  will  be  used  to 
analyze  the  problem  previously  outlined.  It  is  apparent  by  reference  to 
the  boundary  conditions  listed  as  equations  (13)  that  in  equation  (15) 
the  value  of  wx  is  the  known  quantity  and  Avx  is  the  unknown.  Hence, 
equation  (15)  is  an  integral  equation  which  must  be  inverted  in  order 
that  the  solution  can  be  written.  Such  an  inversion  is  not  difficult  if 
the  value  of  wx  is  known  everywhere  in  the  interval  -s1<y1<s1.  In 
the  present  case,  however,  there  is  a  subinterval  ri2  ^  ^  "^l2  in 

which  Wjl  is  not  specified,  and  further,  in  which  Avx  is  not  neces¬ 
sarily  zero  (due  to  the  presence  of  a  trailing  vortex  sheet).  It  will 
be  shown  in  the  subsequent  development  that  the  assumption  that  Avx 
is  zero  in  this  interval  (i.e.,  no  vortices  are  shed  by  the  wing  ahead 
of  the  interval)  yields  a  nearly  constant-chord,  swept-back  wing;  with 
such  a  restriction  the  inversion  can  again  be  performed. 

Given  the  inversion  of  equation  (15),  it  is  possible  to  write  both 
wx  and  vx  for  certain  portions  of  the  real  axis.  All  along  this  axis 
the  functions  wx  and  vx  are,  of  course,  real.  Hence,  if 

fUi)  =  v1(y1,z1)-iw1(y1,z1)  (l6a) 


then  by  analytic  continuation 

f(ll)  =  v(|1,0)-iw1(|1,0)  (l6b) 

Therefore,  the  inversion  of  equation  ( 15) ,  together  with  equation  (l6), 
gives  sufficient  information  to  determine  the  induced  velocities  through¬ 
out  space. 
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Particular  solution  for  the  nearly  constant -chord  wing.-  Adopt  the 
notation  Av ^  equals  Avx  in  the  region  of  the  yx  axis  representing 
the  body  or  fuselage  in  the  £1  plane;  Avia  equals  the  value  of  Av1 
in  the  region  of  the  yx  axis  representing  the  space  between  the  fuse¬ 
lage  and*  the  wing;  and  Avlw  equals  the  value  of  Av±  in  the  region  of 
the  y±  axis  representing  the  wing  plan  form.  Then  if  Avia  =  0  (the 
calculation  of  the  trailing-edge  shape  corresponding  to  such  a  choice 
will  be  presented  later),  equation  (15)  becomes 


w 


l 


2* 


-tx 


Av 


1W 


yi-ya 


-Si 


ri 


-*1 


Avlb 

yi-ya 


dy2- 


l_ 

2« 


_^1S£ 

yi-y2 


dy2 


(17) 


Since  the  airplane  is  laterally  symmetrical,  the  span  loading  is  sym¬ 
metrical  and  A9i(yi)  =  AcPi(-yi).  Therefore,  Avi  has  the  property 
Av1(y1)  =  -  Avi(-yx).  By  means  of  this  relation  for  Avi  and  the  addi¬ 
tional  change  in  notation 


’ll  «  yi 


2  1 


^2 


(18) 


equation  (17)  can  he  written 


=  - 


AvibdTj2 


1_  P1  Avxwdria 

2*  J  2  ni-Tla 
^  1 


(19) 


Equation  (19)  will  now  be  inverted  under  the  condition  that 
Wx  =  -  V0a  for  0<T)i<ri2  and  for  tx2<T]x<Si2,  and  under  the  addi¬ 
tional  condition  that  (Av^w)  2  =  (Avia)  2  =  0,  which  amounts  to 

T^tl 

assuming  the  Kutta  condition  along  the  wing  trailing  edge  (see,  e.g., 
reference  4) .  This  inversion  is  accomplished  by  a  double  application 
of  the  following  solution  (see  appendix  B):  If 


f(ii)  -  - 


1_ 

2it 


b 


a 


Av1(rig)dT]2 


(20a) 


then,  under  the  condition  that  Avx( a)  =  0 


Av1(t)1)  = 


gjjTg)  /b-i lg 

J  r)2~a 


dTia 


(20b) 
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Now  write  equation  (19)  in  the  form 


Si2  Ti2 

_V0a  +  i.  P  m  _  1_  r  AvlbdTlg 

2«J  nx-i2  2*J  %-n2 

4-  n 


and  then,  since  Av1-b(o)  =  0  by  reasons  of  symmetry,  apply 
equation  (20b).  For  0<T)1<r12,  there  results  the  expression 


Avib  =  -  2V0a 


'  %  1 
ri2-%  +  * 


%  rSl  Aviw  /vn 


ri  J  Ti 


Substitute  equation  (21)  back  into  equation  ( 19) ,  reverse  the  order  of 
integration  and,  for  t12<rj1<si2,  there  results 


•V0a  =  -  - 


i_  /  rsi  5(^2) diis 

2n  V  m-ri2  J  _  TI1-TI2 


where 


g(0  =  Av1w(ti2) 


Vri 


Again  apply  equation  (20b),  this  time  to  equation  (22).  In  this  way 
g(rj1)  can  be  shown  to  satisfy  the  relation 

g(%)  =  -  2V0a  /hlhf  I 

V  H  -t|x 

and  equating  this  expression  to  equation  (23)  gives 


AviW  =  -  2VQa 


/  Tli(Tl1-t12) 
(s12-Ti1)(il1-r12) 


,  tj_2  <  <SX2 


(25a) 
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The  results  given  by  equations  (25)  can  be  extended  to  other  points 
in  space  by  analytic  continuation.  Equations  (l6)  indicate  the  necessary 
procedure.  Hence,  since  (v^  _  =  Av-l/2, 


3When  the  method  is  applied  to  a  value  of  wx  which  has  some  given 
variation  with  r\1  there  results 
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Equation  (26)  has  several  branch  points  so  it  is  not  uniquely  defined 


without  specifying  the  cut  from 
the  £3,  plane  (see  sketch  (c)). 

In  the  upper  half  of  the  li  plane 
varies  between  0  and  rt  and 
in  the  lower  half,  between  0 
and  -rt.  Notice  that  when  ii  is 
at  a  very  large  distance  from  the 
origin  in  any  direction  the  mag¬ 
nitude  of  the  term  on  the  right- 
hand  side  of  equation  (26)  tends 
to  zero,  so  that  the  boundary 
conditions  at  infinity  are  satis¬ 
fied.  It  is  evident  that  the 
other  boundary  conditions  in 
equations  (13)  are  also  satisfied. 


to  si  along  the  real  axis  in 


(C) 


The  Trailing  Edge 


Special  trailing-edge  shape.-  Equation  (26)  is  a  solution  to 
Laplace's  equation  and  represents  the  flow  around  a  wing  and  body. 
However,  the  plan  form  of  the  wing  has  not  as  yet  been  evaluated, 
although  it  has  been  fixed  as  that  which  makes  the  value  of  Avi  vanish 
in  the  region  between  the  wing  and  the  body.  Since  Avi  is  the 
gradient  of  Acpx  in  the  yi  direction,  and  further,  since  (ATi)rjn  ^ 

(the  value  of  ACpi  at  the  wing  trailing  edge)  equals  the  total  circu¬ 
lation  F1  about  a  given  chordwise  section,  this  amounts  to  the  same 
thing  as  assuming  that  there  are  no  trailing  vorticies  between  the  wing 
and  the  body.  It  is  a  further  consequence  of  such  an  assumption  that 
the  span  loading  ahead  of  this  region  is  a  constant  for  r02<y02<t02. 
The  configuration  which  will  produce  such  a  flow  must  now  be  determined. 
In  particular,  if  the  leading  edge  is  taken  to  be  a  straight  line,  the 
equation  for  the  trailing  edge  is  unique  and  needs  to  be  expressed. 

One  of  the  simplest  ways  of  finding  the  shape  of  the  trailing  edge 
is  to  find  ATt.E.  from  equation  (26)  and  solve  for  t  as  a  function 
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of  s  for  a  fixed  value  of  rQ  and  AT^.e.*  The  constant  represent¬ 
ing  is  the  value  of  the  jump  in  potential  at  the  point  P  in 

sketch  (d).  Here  Acp  is  known  (see,  e.g.,  reference  3)  since  there  is 

no  gap  to  make  its  solution  indeter¬ 
minate. 


Consider  an  arbitrary  section,  as 
AA  in  sketch  ( d) ,  downstream  of  the 
point  P.  The  value  of  Aviw  at  such 
a  section  is  given  by  equation  (25a) 
and  the  solution  for  Acpiw  follows  by 
definition  and  is 


AT1W=  ~2V„a 


y  2-t  1 
^2  1 


(sia-y22)(y/-ri2) 


(27) 


Equation  (27)  is  an  elliptic  integral 
which  can  be  easily  reduced  by  means  of 
the  substitutions 


ki 


2 


Si2-t*a  k . 

si2'r12  ’  1 


(28) 


and  by  using  the  Jacobian  elliptic  functions  defined,  in  this  case,  by 
(si2-ti2)sn2u  =  si2-y!2,  cn2u  =  l-sn2u,  dn2u  =  l-kx2sn2u 
to  the  form 


where 


^(^1,^  1)  -kjj  F(kx,  +x) 


(29) 


and  where  the  incomplete  elliptic  integrals  E  and  F  are  defined  in 
the  list  of  symbols  (appendix  A).  Equation  (29)  reduces  to  the  results 
given  in  references  3  and  k  when  there  is  no  gap  or  no  body,  respectively 
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At  the  trailing  edge  y ^  =  tj_,  and  equation  (29)  becomes 

(ACP1)t.e.  =  2V0a  ( 

where  the  elliptic  integrals  are  now  complete.  Transform  this  to  the 
physical  plane,  using  equation  (6),  and  set 


t(sW) 


9  1  “  J  l“k( 


then  there  results 


A^T.E.  =  2V0a  (i~-)  (E0-k0’2Ko)  (33a) 

At  the  juncture  of  the  fuselage  and  the  wing  trailing  edge  (the  point  P 
in  sketch  (d)),  s  equals  mc0  and  t  equals  rQ,  so  that 
equation  (33a)  reduces  to 

(Acp)T.E.  =  2V0amco  1-  (33b) 

As  was  pointed  out,  the  solution  for  the  equation  of  the  trailing  edge 
can  be  obtained  by  equating  these  two  values  of  (a<P)ij>  £  .  Hence 


Ik 
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and  since  s  equals  mx 

,  .  _  ( “* 

°  \  m2x2-r02/  t 

From  equation  ( 34)  the  solution  that  gives  the  correct  trailing-edge 
shape  can  he  written 


t  = 


(35) 


and  from  the  definition  of  ko? 


mx 


t£-r02  +  J (tg-rQ2)2+4kQ ' 2t2rQ2 

2ko't 


(36) 


If  k0T  and  r0/mc0  are  fixed,  t/mcQ  is  determined  from  equation  (35) 
and  a  fixed  k0f  ,  r0/mco  and  t/mc0  determines  x/cQ  from 
equation  (36).  Hence,  it  is  relatively  easy  to  calculate  numerically 
the  shape  of  the  trailing  edge. 

Sketch  (e)  shows  the  shape  and  position  of  the  trailing  edge  when 
the  wing  leading  edge  is  swept  back  45°  and  the  radius  of  the  fuselage 


is  31*6  percent  of  the  extended  root 


chord,  cQ*  (A  dimensionless 

coordinate  system  is  chosen, 
however,  so  that  the  results 
can  be  used  for  various  values 
of  m  and  cQ.)  Shown  also, 
for  comparative  purposes,  is 
the  position  of  the  trailing 
edge  when  there  is  no  fuselage 
the  condition  in  both  cases 
being,  of  course,  that  no 
vortices  trail  back  in  the 
region  directly  behind  the  wing 
Table  1  presents  coordinates  of 
the  trailing  edge  for  several 
values  of  r0/mc0. 
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Sketch  (f)  indicates  the 
variation  of  the  local  chord 
along  the  span  for  r0/mc0 
equal  to  0.6,  0.316,  and  0. 

It  is  apparent  that  the 
effect  of  the  body  is  to  make 
the  trailing  edge  more  nearly 
that  of  a  constant -chord  wing. 
The  asymptotic  value  of  the 
wing  chord  is  given  by  the 
equation 


where 

lim  f  t'N 

C°°  “  00  \x"  m/ 

More  complete  results  are 
given  in  table  2. 


Other  trailing-edge  shapes.-  The  procedure  just  presented  can  be 
generalized  and  used  to  calculate  trailing-edge  shapes  corresponding  to 


arbitrary  span-loading  curves.  Suppose  that  the  span  loading  in  the 
transformed  (or  £x)  plane  is  represented  by  a  power  series  in  y1  in 
the  interval  ri2<yi2<  tx2.  Then  the  circulation  in  this  interval  can 
be  written 
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1_ 

2n 

The  left-hand  side  of  the  latter  equation  varies  with  tji  in  a  given 
manner  depending  on  the  b^s  in  the  expression  for  T.  Hence,  the 
equation  can  be  considered  as  identical  to  equation  (19),  the  left  side 
being  regarded  as  an  effective  wx  in  equation  (19).  The  analysis 
succeeding  equation  (19)  can  now  be  repeated  in  terms  of  the  equivalent 
wp.  There  results  an  equation  for  the  trailing  edge  which  depends  on 
the  bn*s. 

By  the  process  outlined,  both  the  trailing-edge  shape  and  the  span¬ 
loading  curve  have  been  expressed  in  terms  of  m  +  1  constants.  By 
varying  the  number  and  magnitude  of  these  constants,  a  large  class  of 
trailing-edge  shapes  can  be  obtained. 


AvibdT); 


1_ 

2n 


Avlwa-n 


T\  -T) 
*1  '2 


*1 


Wi 


2x  u 


rr 


AviadTl2 

T]  -'H 
*1  *2 


The  Wing  Area 


Having  found  the  shape  of  the  trailing  edge  by  the  methods  outlined 
in  the  preceding  section,  it  is  now  possible  to  determine  the  area  of 
the  wing.  Denote  this  area,  region  1  in  sketch  (g),  as  and  the  area 
of  region  2,  shown  also  in  the  sketch,  as  S2.  It  is  evident  that  the 
sum  of  these  two  areas  is  simply 


(9) 
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If dimensionless  system  based  on  the  length  c0  is  adopted,  one  can 
write  for  the  total  area  (i.e.,  both  panels,  see  the  shaded  area  in 
sketch  (hi))  of  the  wing  the  equation 


where  t/mc0  is  given  numerically  as  a  function  of  x/cQ  in  table  1 
(y/mc0  in  the  table  representing  t/mc0) •  Numerical  relations  between 
the  parameters  S/mc0  ,  r0/mc0,  s0/mc0,  and  t0/s0  are  presented  in 
figure  1. 

The  area  of  a  wing  with  another 
kind  of  tip  shape  can  be  readily 
evaluated  once  the  particular  tip 
shape  is  specified.  For  example,  the 
area  of  the  wing  shown  in  sketch  (h2) 
can  be  calculated  by  subtracting  a 
rectangular  area  (given  by  the  sum  of 
the  two  triangular  regions  labeled  3 
in  sketch  (h2))  from  the  area  of  a 
sketch  (hi). 


Downwash  Behind  Wing 

The  equation  for  the  downwash 
behind  the  wing  and  in  the  z  =  0 
plane  follows  immediately  from 
equation  (26).  In  the  transformed 
l-L  plane  the  value  of  wx  is 

W1  =  -  V0ct,  0<yi2<r12;t12<y12<si2 

(38a) 

and 


(38t>) 
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In  order  to  transform  this  value  to  the  physical  plane,  care  must  be 
taken  to  go  backwards  through  the  boundary  conditions  in  the  proper 
order.  Equations  (38)  represent  the  solution  for  the  boundary  conditions 
presented  in  equation  (13).  To  find  the  solution  for  the  conditions 
given  by  equation  (l2)  a  free  stream  V0<x  must  be  added.  Thus,  in 
mathematical  notation. 


(wi)  +  V0a 

x  13 


where  the  subscripts  12  and  13  refer  to  the  boundary  conditions  satis¬ 
fied.  Finally,  to  find  the  downwash  in  the  physical  plane,  the  transfor¬ 
mations  given  as  equations  (9)  and  (5)  must  be  employed  and  the  free 
stream  subtracted  so  that 


w(y)  = 


(wi)i3  +  v°a 


y2-^ 


-  V„a 


where  (w-i)  now  becomes 

1  13 


(wx)13  =  -  VDa,  0<  y2<  r2;t2<  y2<  s2 


and 


(wl)13  =  -  VQa 


1  - 


y2-r  2 


■r0 


(t 


y*+r02  /B\  /(y*t2-r0«)(t2-y*) 


(B2y2-r04)(sV) 


,  r2<  y2<  t2 


Combining  and  simplifying,  one  finds 

w(y)  =  -  Voa,0<y2<r2jt2<y2<  s2  (39a) 

and 


w(y) 


VGa 


2  2 
y  +r0 


/  2,2  4 \ /  2  2 \ 

(y  t  -r0  )(t  -y  ) 
(s2y2-r04)(s2-y2)  J 


,  r2  <  y2<  t2 


(39b) 
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The  accompanying  sketch  shows  the  variation  of  -w/V0a  in  the 
intervals  for  which  it  has  heen  given.  If  no  wing  is  attached  to  the 


body  (or  if  the  gap  is  very  large)  the  fluid  at  the  side  of  the  body  is 
moving  upward  at  a  speed  equal  to  that  at  which  the  body  is  moving  down¬ 
ward.  The  presence  of  the  wing  restricts  this  motion  and  as  the  wing 
panel  approaches  the  body  the  air  in  the  gap  is  forced  more  and  more  to 
move  downward  with  the  wing  and  body.  The  dotted  lines  in  the  sketch 
show  the  variation  of  -w/V0a  if  no  body  is  present,  that  is,  if  rQ 
equals  zero. 


Chordwise  Load  Distribution 


Loading  on  the  wing.-  The  loading  on  the  wing  can  be  calculated 
by  means  of  the  linearized  equation  for  the  loading  coefficient.  This 
equation  can  be  written  ♦ 


2Au  _  2 

v7  '  ~ 


(4-Oa) 
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It  is  somewhat  easier  to  calculate  the  loading  if  the  derivative  of  Acpw 
is  taken  in  the  plane.  If  A<P1V  is  considered  to  he  a  function 

of  the  two  independent  variables  yx  and  sL,  equation  (40a)  can  he 
modified  slightly  to  read 

/^p  \  2_  ggPiy  dfli  is 
\  q  V0  dsx  ds  dx 


(40h) 


The  value  of 
thus  dsi 


^Piw  ___  be  stained  hy  differentiating  equation  (29), 


can 


Ba<P 


iw 


dsx 


=  2V0a 


J~  2  "  2 


E(ki,+i)-k1»2F(k1,ti)  + 


bx  -ra 

^(ki,^) 


] 


L 


-  k- 


aF(kl’*d  -  zvfo^x) 


ds 


ds 


which  hecomes 


dA<Piw 


els 


-  =  2Voa 


ysi2-ri 


/-— -%-g) 

N  /  l-ki2^2/ 


2  /  l-ti2  gj[x 
dsx  1  7  l-kj2^2  dsi 


dk. 

ki  — i  +k 


(*U) 


The  terms  dk-j/ds^L  and  d+i/dsi  both  involve  dt-^ds-L  which  is 
proportional  to  the  slope  of  the  trailing  edge  in  the  transformed  £1 
plane.  This  latter  derivative  can  be  readily  obtained  from  equation  (^l) 
itself  since  the  value  of  (Ap/q)v  and  hence  dAcpiv/dsi  must  be  zero 
on  the  trailing  edge,  that  is,  where  yj_  equals  and  ^  equals  1* 

This  yields  the  relation 

dtl  =  siEi 
cL^l  tj^Kj^ 


(42) 
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by  means  of  which  the  identity 


— i  =  S1  /tA  1-+!2 

dsi  s xs-rxz  \  ki2^  +  +~k~2 


dsi  si2-ri2 


Ei-ki*  KX 


there^esultg-for^th^loading0113  ^  (“3)  ^  'qU8tl0n  <to><  ■“* 


1  /Ap 


m  V  qa 


‘*7^;.  ['(‘-♦J-g'lMJ  ♦ 


^  yu-Oci-ii2**)  _ 

dsi/ds^U2  r^Ws2  ^  ?  plane>  one  findsj  finally,  si 


since 


-  (^-p)  -  ^ 
m  V  qa-4 


2  2 
s  +rQ 


E(ko,Hr0)  -  -  F(ko,+0)  + 


ffjjfgof)  /(y2- t2)(ygt2-rQ4) 
yt( S2-r02)  J  ( S2_y2) ( s2y2_rQ4) 


where  ko  is  defined  by  equation  (33)  and 


*Q  =  !  /(s2-y2)(s2y2-r04) 

y  /(s£-t2)(82ta-r0‘)  (45) 

_  In  the  special  cases  when  there  is  no  body  or  when  the  wing  is 

ref M Uat^°?J44  agreeS  W±th  the  results  Presented  in 

* ®S.^].and  (3)>  respectively.  A  discussion  of  the  chordwise 
load  distribution  over  the  wing  will  be  given  at  the  end  of  this  section 


22 


NACA  TN  2554 


Loading  on  the  “body.-  The  variation  of  the  load  distribution  over 
the  body  can  be  calculated  in  much  the  same  way  as  that  over  the  wing. 
It  is  first  necessary,  therefore,  to  find  the  jump  in  potential  between 
points  directly  opposed  above  and  below  the  z  =  0  plane.  In  the 
plane  this  difference  follows  immediately  from  equation  (25b)  just  as 
equation  (27)  was  written  for  A<Piw  Hence, 


Acpib  =  -  2V0a  /  y2 


y22-ti2 


(s!2-y22)(y22-ri2) 


dy2 


r1 

2V0a  /  y2 

r 


,2  2 
ti  -y2 


(si2-y22)(ri2-y22) 


<iy2 


(46) 


The  first  of  the  integrals  in  equation  (46)  has  already  been  evaluated 
and  the  second  can  be  reduced  by  means  of  the  transformation 

(ti2-y22)  sn2^  =  rx2-y22 
After  some  manipulation,  equation  (46)  becomes 


A<Pib  =  2V0a 


2V0a 


E1-E:(ki,t3)»ki'2F(kl,*3)-k1'2K1 


y 


si2-ri2  + 


(si2-yi2)(ri2-yi2) 
ti2-yi2 


where  kx  is  defined  by  equation  ( 28)  and 


(47) 


**  = 


I  2  „  2 
ri  -yi 

4.  2  v  2 

ti  -yi 


(48) 


Using  the  equations  (42)  and  (47),  one  can  write  for  the  loading  on  the 
body  (after  differentiation  and  simplification  of  equation  (47)) 
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1  (4B  )  _  4 

m  \qa  J 


-2  F(k0,*2)-E(ko,t2)  + 
Ko 


2ko2( s2-r02)t 


ro2-y2 


V  C(t2+r02)2-4y2t2] [ ( s2+r02)2-4y2s2] 
where  kQ  is  defined  hy  equation  ( 33)  and 


+2  =  2t 


ro2-y2 

( t2+r02) 2-4y2t2 


Again  equation  (49)  agrees  with 
previously  known  results  in  the 
limiting  cases  when  the  body 
vanishes  or  when  the  wing 
becomes  triangular. 

Discussion  of  the  chord- 
wise  loading . -  Equations  (44) 
and  (49)  form  the  basic  results 
of  part  I  of  this  report. 

Graphs  of  the  loading  coeffi¬ 
cient  for  a  wing  alone  and  for 
a  wing -body  combination  (r0/mc0 
equals  O.316)  are  shown  in 
sketch  (jl).  The  results  for 
the  case  of  zero  body  radius 
could  have  been  obtained 
directly  from  reference  k . 

They  are  shown  here  for  the 
purpose  of  a  qualitative  com¬ 
parison.  Unfortunately,  the 
load  distributions  on  the  two 
wings  cannot  be  compared  quan- 
tatively  on  the  basis  of  equiv¬ 
alent  plan  forms  since  the 
trailing-edge  shapes  differ 
significantly.  The  variation 
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of  the  loading  along  the  center  line  of  the  body  is  shown  in  sketch  (j2). 


On  the  basis  of  the  load  distributions  presented  in  references  3 
and  4,  the  qualitative  variation  of  loading  shown  in  sketch  (jl)  is 
obvious.  That  is,  the  loading  falls  steadily  from  its  infinite  peak  at 
the’  leading  edge  to  zero  at  the  trailing  edge.  On  sections  which  are 
cut  by  the  Mach  wave  from  the  trailing-edge  fuselage  juncture,  the  slope 
of  the  curve  is  discontinuous. 

The  change  in  the  load  distribution  brought  about  by  the  presence 
of  a  wing  tip  is  the  same  for  a  wing-body  combination  as  for  a  wing 
alone.  The  behavior  of  the  loading  in  the  vicinity  of  a  tip  has  a 
straightforward  explanation  in  terms  of  the  trailing  vortex  sheet. 

Thus,  if  the  wing  is  cut  off  along  a  line  perpendicular  to  the  free- 
stream  directi cn,  the  vortices  which  were  bound  in  the  wing  all  turn  and 
trail  backwards  with  the  same  distribution  in  strength4  as  they  had  when 


4This  assumes,  of  course,  that  the  vortices  have  not  begun  to  roll  up  to 
any  significant  extent. _ 
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crossing  the  last  spanwise  section  of  the  wing 
(see  sketch  (k)).  Since  the  vertical  induced 
velocity  along  the  last  spanwise  section  was  made 
constant  (by  finding  the  appropriate  solution  to 
the  integral  equation) ,  it  must  also  be  constant  - 
and,  in  fact,  the  same  constant  -  everywhere  in 
the  vortex  wake.  Hence,  if  a  flat  surface  having 
the  same  angle  of  attack  as  the  wing  is  inserted 
anywhere  in  the  wake  it  will  in  no  way  disturb  the 
flow  and  consequently  there  will  be  no  loading  on 
such  ^a  surface  ( just  as  there  is  no  loading  on  the 
vortex  wake  itself) .  The  loading  is  zero,  therefore, 
on  the  tip  regions  marked  2  in  sketch  (k);  the  load¬ 
ing  in  the  regions  marked  1  being  given,  of  course, 
by  equation  (44). 

It  is  interesting  to  see  how  the  distribution 
and  magnitude  of  the  loading  given -by  this 
(slender  wing)  theory  compare  with  linearized 
theory  results  at  some  Mach  number  other  than  1. 

The  differences  caused  by  considering  Mach  numbers 
other  than  1  depend,  of  course,  on  whether  or  not 
the  new  Mach  number  is  subsonic  or  supersonic. 

This  discussion  must  be  limited  to  a  comparison 
with  supersonic  Mach  numbers  only,  since  theo¬ 
retical  chordwise  load  distributions  over  swept- 
back  wings  flying  at  high  subsonic  speeds  are  not 
available .  The  change  in  the  loading  brought  about  by 
increasing  the  speed  can  be  divided  into  two  parts: 
one,  a  change  caused  by  the  rotation  of  the  Mach 
lines  which  form  the  boundaries  of  the  various 
regions  in  each  of  which  the  shape  of  the  loading 
curve  takes  widely  different  forms;  and  the  other, 
a  change  in  the  magnitude  of  the  loading  within 
each  of  these  regions. 


Sketch  (l)  indicates  these  effects.  Thus,  on  the  wing  flying  at 
supersonic  speeds  the  sharp  drop  in  loading  occurring  at  a  critical  Mach 
line  moves  farther  back  along  the  chord  from  point  b  to  point  a  in 
sections  AA  and  BB  shown  in  the  sketch.  This  causes  a  considerably 
higher  value  of  the  loading  for  the  supersonic  wing5  in  regions  1  and  2. 

A  similar  effect  occurs  on  the  body  traveling  at  a  supersonic  speed  where 
now,  however,  the  traces  of  the  Mach  lines  are  no  longer  straight  but, 


5 Solutions  showing  the  effect  of  crossing  critical  Mach  lines  on  a 
swept-back  supersonic  wing  are  given  in  references  6  and  7. 
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due  to  the  curvature  of  the  body,  form 
helices.  Region  3  in  the  supersonic  case 
would  he  a  region  of  zero  loading  and  region  4 
would  he  a  region  of  high  loading  relative 
to  the  sonic  value. 


The  relative  magnitude  of  Ap/qa  within 
the  various  areas  hounded  hy  the  wing  edges 
and  the  pertinent  Mach  lines  changes  as  the 
reference  section  moves  outboard  along  the 
wing  span.  Along  inboard  sections  ahead  of 
region  1  in  the  sketch  (i.e.,  ahead  of  the 
sonic  Mach  line  from  the  trailing  edge  root) 
the  loading  on  the  sonic  wing  is  higher  than 
that  on  the  supersonic  wing.  It  is  well 
known,  for  example,  that  in  the  case  of  a 
triangular  wing  without  body,  slender -wing 
theory  gives  a  loading  E  times  the  loading 
obtained  at  a  supersonic  Mach  number  (where 
E  is  the  complete  elliptic  integral  of  the 


second  kind  with  modulus 


V^l-m2p2 


and  is  given 


closely  by  1  +  i  m2P2  (in  —  -  -  ^  for 

2  \  mp  2  J 

small  values  of  mp) .  On  the  other  hand, 
along  sections  farther  outboard,  the  magni¬ 
tude  of  Ap/qa  on  the  sonic  wing  must  become 
lower  than  that  on  the  supersonic  wing.  This 
follows  immediately  from  simple  sweep  theory,  * 
since  the  component  of  velocity  normal  to  the 
leading  edge  is  closer  to  the  speed  of  sound 
than  that  for  the  sonic  wing.  In  fact,  it  is 
easy  to  show  that  at  distances  far  enough 
outboard  so  that  simple  sweep  theory  applies6 
the  supersonic  wing  has  a  loading 


,2 io2  a-1/2 


1  -  wr$*J  '  times  that  obtained  from 
slender  wing  theory. 


By  an  application  of  the  above  consider¬ 
ations,  it  is  possible  to  obtain  an  estimate 
of  the  absolute  value  of  the  loading  on  a 
wing-body  combination  at  supersonic  Mach 


6 Sketch  (jl)  indicates  the  manner  in  which  the  loading  approaches  that 
given  by  simple  sweep  theory  as  the  reference  station  moves  outboard. 


NACA  TN  2554 


27 


numbers.  Another  manner  in  which  the  results  of  slender  wing  theory  can 
he  extended  to  Mach  numbers  other  than  1  (or  to  plan  forms  which  are  not 
sufficiently  slender)  is  to  form  the  ratio  of  the  resulting  values  for 
the  wing  plus  body  to  those  for  the  wing  alone  and  apply  this  ratio  to 
solutions  for  the  same  wing  or  body  at  the  required  flight  Mach  numbers 
( or  slenderness  factor).  As  has  already  been  mentioned,  the'  formation 
of  such  a  ratio  for  the  load  distribution  is  not  possible  from  the  solu¬ 
tions  presented  herein  since  the  wing  trailing  edges  change  to  a  certain 
extent  with  the  addition  of  the  body.  It  is  reasonable  to  expect,  how¬ 
ever,  that  a  ratio  of  the  integrated  loading  characteristics  (i.e.,  lift, 
drag,  and  pitching  moment)  formed  by  dividing  the  result  for  a  wing-body 
combination  by  those  for  a  wing  alone  will  be  useful  in  estimating  the 
interference  effects  even  if  the. wing  trailing  edges  differ  slightly. 


Aerodynamic  Characteristics 


The  results  developed  in  the  preceding  section  can  now  be  converted 
into  forms  which  represent  the  aerodynamic  characteristics  of  the  wing 
and  body.  Hence,  the  following  will  present  the  span  loading,  average 
chord  loading,  lift,  drag,  pitching  moment,  and  center  of  pressure  for 
the  wing -body  combination. 

Span  loading.-  The  development  of  the  span  loading  on  the  wing  and 
body  will  be  considered  separately.  First  the  span  loading  on  the  wing 
can  easily  be  determined  from  the  value  of  A<P  given  in  a  preceding 
section.  Thus 

T-  If 

wing  plan  form 

and  since  ACp  at  the  leading  edge  is  zero 

^  =  %  f  (a<p)t.e.  dy 

span 

where  (A9)ip>E>  is  the  value  of  A9  on  the  trailing  edge  of  the  wing. 
Since  (A<p)T>E>  also  represents  the  total  circulation  about  the  wing 
chord,  there  results  for  the  circulation  Tw  developed  by  the  wing  and 
the  total  wing  lift  Ly- 
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rw  =  (A^T.E. 


The  equation  for  Tw  can  now  be 
determined  from  equation  (29). 

Between  y  -  r0  and  y  =  tQ  (see 
sketch  (m))  the  value  of  Tw  is  a 
constant  (this  being  the  condition 
by  which  the  shape  of  the  trailing 
edge  was  determined) .  Between  y  =  tQ 
and  y  =  s0,  is  given  by  the  value 
of  AT  along  the  section  AA  since 
there  is  no  loading  between  this 
section  and  the  trailing  edge.  Hence 


rw  =  2V0ac0m 


(52a) 


and 


E(ko,Vo)-ko'  F(k<>,*o) 


t0<  y<  s0 


(52b) 


where  k0  and  are  defined  in  the  table  of  symbols. 

Some  care  must  be  taken  in  order  to  find  the  span  loading  on  the 
body.  Since  we  are  concerned  here  with  the  loading  developed  behind  the 
wing-leading-edge  fuselage  juncture,  it  is  necessary  to  subtract  the 
value  of  ACj>b  at  this  station,  shown  as  station  AA  in  sketch  (n),  from 
ATft  at  station  BB  also  shown  in  sketch  (n) .  For  the  total  span  loading, 
then,  it  will  be  necessary  to  add  to  this  value  the  load  accumulated 
on  the  nose  of  the  body.  Denote  by  (i TD)0  the  increment  of  circulation 
developed  by  the  nose  of  the  body  and  by  (Fb)i  the  increment  of  circu¬ 
lation  developed  behind  the  wing-leading-edge  fuselage  juncture,  that 
is,  between  stations  AA  and  BB.7  Hence> 


7Slender  wing  theory  gives  zero  loading  behind  station  BB  as  long  as 
the  trailing  vortex  pattern  does  not  vary. 
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(rb)1  =  (Acp)T.E.  -  (Acp)L.E. 

»  (53) 

(•Lb)  1  =  PoVq  ['  (A)^  cLy 


The  value  of  (Acp)^^.^  see  sketch  (n), 
can  be  obtained  from  equation  (47)  by 
setting  tx  and  si  equal  to  n* 

By  transformation  of  the  result  into 
the  physical  plane ,  one  obtains 


The  span  loading  on  the  body  is  then 
given  by 


where  kQ  and  are  given  in  the  table  of  symbols. 
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r 


VoaCofT) 


table  2.  Typical  results  are  shown 
and  0.316  in  sketch  (p). 


Sketch  (o)  shows  the 
variation  of  the  span  loading 
over  the  wing  and  body  for  a 
body  radius  factor,  r0/mc0. 
equal  to  0.316  and  a  wing 
semispan  factor,  s0/mc0 
equal  to  1.7. 

Section  lift.-  The  wing- 
section  lift  coefficient  can 
be  calculated  readily  by 
dividing  the  section  chord 
into  the  value  of  the  span 
loading  at  the  same  span  sta¬ 
tion.  Along  sections  not 
influenced  by  the  tip  cut-off 
this  is  especially  simple 
since  the  span  loading  is 
constant.  The  value  of  the 
section  lift  curve  follows 
immediately,  therefore,  from 
a  body  radius  factor  equal  to  0 


0 


tp) 


/ 


2  y/mc0  3 
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Section  drag.-  The  value  of  the  section  drag  can  he  written 

cd  =  acz  +  (cd)  (55) 

where  (c^)s  represents  the  suction  force  at  the  section  leading  edge. 
The  magnitude  of  (ca)s  can  he  evaluated  (see,  e.g.,  reference  8)  by  the 
equation 

(cd)s  (56) 

qc  dy 

where  dF/dy  is  the  suction  force  in  the  free -stream  direction  per  unit 
length  normal  to  the  free  stream,  and  c  is  the  local  chord. 

Define  a  new  set  of  coordinates,  as  shown  in  sketch  (q),  such  that 
yn  lies  along  and  xn  lies  perpendicular  to  the  leading  edge  of  one 
wing  panel.  Then  if 

g  =  ii”  en  ."A* siz  y 

Xn  O  n  Vn  n  ' 


/'/  \  Leading 
/  \  edge7 


the  suction-force  component  F  (positive 
in  the  positive  xn  direction)  in  the 
free-stream  direction  is  given  hy  the 
equation 

f  -  -  G2U,7)  (58) 

dy  Pn 


Now  hy  differentiating  equation  (29)  with  respect  to  x  and  dividing 
hy  2  (to  convert  Au  into  u) ,  there  results 


u  =  mV0a 


,2,.  2 


E(k0,*0)-  -  F(ko,*o)+ 
Kn 


s5(y2-r02) 

yt(s2-rQ2) 


Ay2-t2)(y2t2-r04) 
( s2-y2) ( s2y2-r04) 


y/mcQ  >  1 
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and  similarily  v  becomes 


v , .  foaiv  (1)  / ,  y/MSo> ! 

y2  \t/  v  ( s2y2-r04)  ( s2-y2) 

Since  the  normal  component  un  is  given  by  the  equation 


un  = 


mu-v 


J  1+m2 


there  results 


un  =  m 


J 


l+m‘ 


V0a  f s2+r02 


2  \  S' 


E, 


E(k0,t0)-  F(ko,i0) 

*o 


Vpa  /(y2-t2)  (y^2-^4) 
J  yt/1^  y  (s2-y2)(s2y2-r04) 


m£(s2+r02)(y2-r0_f)_  s(y2+rQ2) 


s2-r02 


,  y/mcQ  >  1 


By  means  of  equation  (57) >  G  can  now  be  calculated.  Hence,  since 
Pn  =  1  /  <\/  1+ni2  y 


G  = 


a(y2+r02)  /(y2-t2) (y2t2-r04) 


yt(l+m2) 1^4  J  2y(y*-r04) 
Finally,  therefore,  the  suction  force  can  he  written 

a2(y2+r02) 


dF 

dy  "  2m  y*t2 


(y^i^Hy^-ro4)  1 


2y(y2-r02) 


(59) 


and,  by  using  equations  (55)  and  (56)>  the  section  drag  coefficient  can 
he  written 
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cd  cl  «(y/mc0)  L1+ 


(?)'][>■(!)•][*•  (?)'(¥) 


a2m  am 


(c/co) 


-  (—1 
\y J 


i-\ — )  (^)l 

y/mc0>l 


(60a) 


In  the  region  where  the  leading-edge  suction  force  cannot  he  affected  hy 
the  trailing -edge  shape,  equation  (60)  reduces  to  the  simpler  form 


cd  ...  CZ 
am 


4 

jt(y/mc0)  ±  r_£\ 

(e/c0)  1- 


ro/mc0<y/mc0<l 


(60b ) 


The  variation  of  the  section  drag  coefficient  is  shown  in  sketch  (r) 
for  two  wings;  one  without  a  "body,  and  the  other  with  a  hody  radius 
factor  equal  to  0*316* 
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Total  lift.-  The  lift  on  the  wing,  the  "body,  and  the  combination 
can  now  he  evaluated  hy  means  of  equations  (51)  and  (53)  since  the 
expressions  for  r  are  given  hy  equations  (52)  and  (54).  The  integra¬ 
tion  required  is  somewhat  involved  algebraically,  hut  the  final  result 
can  again  he  expressed  in  terms  of  elliptic  integrals.  Thus,  defining 
A Q(k,t)  hy 

A0(k,t)  =  ?  KE*(k,f)  +EF'(Mr)  -  KF*  (k,ijf )  (6l) 

(tabular  values  for  A0  can  be  found  in  reference  9),  the  total  lift 
carried  hy  the  wing  is 


4w  _  2  so2~ro2 
qa  s0 


t02+r02 


E2-ko  K2j  -^r, 


i  E0-k0' 
u 


(62) 

Equation  (62)  agrees  with  the  results  presented  in  references  3  and  4 
when  tQ  =  r0  (the  case  of  a  triangular  wing  on  a  body)  and  rQ  =  0 
(the  case  of  no  body),  respectively.  When  sQ  equals  tQ,  that  is, 
when  there  is  no  wing,  L*,  equals  0. 

The  lift  on  the  body  will  he  computed  in  two  parts  just  as  was  the 
span  loading  on  the  wing:  the  lift  on  the  portion  of  the  body  behind  the 
wing-leading-edge  fuselage  juncture  (Lb^,  and  the  lift  on  the  nose  of 
the  body  (l*,)0.  It  is  a  well-known  result  of*  Munk r s  airship  theory  that 
the  lift  on  the  pointed  nose  is  just 


(Lb)0 


2rtrQ2 


(63) 


and  is  independent  of  the  shape  of  the  nose.8  The  value  of  the  lift  in 
the  vicinity  of  the  wing  follows  from  the  integration  of  equation  (54) 
according  to  equation  (53).  The  total  lift  on  the  body  can  then  be 
written 


8 In  this  report,  it  is  assumed  that  the  nose  is  always  ahead  of  the  wing, 
that  is,  the  portion  of  the  body  on  which  the  wing  is  mounted  is  every- 
where  a  circular  cylinder. 
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Lb  (l^)0+(Lb)i 


qa 


qa 


Sq2-^o2 

so 


krQ  { Eo-ko*  K( 


•f  2 

tn  +r  n 


E2-ko  K2 


7t 


Ao(k2^4)  -  2*r02 


(64) 


Setting  tD  =  rG,  one  finds  the  result  given  in  reference  3  for  a  trian¬ 
gular  wing  mounted  on  the  cylindrical  portion  of  a  pointed  body. 

If  tc  =  rQ  =  sG,  the  wing  disappears  and  equation  (64)  reduces9  to 
equation  (63).  If  rQ  =  0,  reduces  to  zero. 

Finally,  the  sum  of  equations  ( 62 )  and  (64)  gives  for  the  total 
lift  of  the  wing-body  combination,  including  the  nose  of  the  body. 


—  =  2«(^±r°t  -  +  a) 

qa  \  Sq2  t02  / 

(65) 

Sketches  (s)  and  (t)  show 
the  total  lift  on  various  wing- 
tody  combinations  together  with 
its  division  into  the  component 
parts  carried  separately  by  the 
wing  and  body.  Various  lift 
coefficients  ,  depending  on  the 
choice  of  the  reference  area, 
can  be  formulated  by  means  of 
the  area-span  relationship 
given  in  figure  1. 

Total  drag.-  In  general, 
the  vortex  drag  can  be  calcu¬ 
lated  by  finding  the  momentum 
transport  through  a  plane  perpen¬ 
dicular  to  the  x  axis  and 
located  infinitely  far  behind 
the  airplane.  In  slender  wing 
theory  the  calculation  of  the 
total  drag  is  simplified  in  two 


9  Note  that  AQ(k,l)  =  1. 
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ways:  first,  the  vortex  drag  becomes  the 
total  drag  (neglecting,  of  course,  vis¬ 
cosity),  and  second,  in  the  calculation  of 
this  drag  the  reference  plane  can  be  located 
immediately  behind  the  airplane  since  the 
flow  there  is  the  same  as  it  is  infinitely 
far  back. 

Hence,  a  momentum  balance  gives  for 
the  drag 


D  =  -  pQ  P  w0A<Pdy+poro  /  (vr)  (<P) 


d0 


r=r, 


r=r, 


(66a) 


where  w  is  the  value  of  the  vertical 
induced  velocity  behind  the  wing  in  the 
z  =  0  plane.  It  is  more  cqnvenient  to 
perform  this  integration  in  the  plane. 

Equation  (66a)  can  be  put  in  the  form 


% 


D 


yi 


-Ai 


v  Atp^dy i~P0 


rTl 

/" 


Acpidyj. 


♦ 

(66b) 


•where  w  and  A<Pi  axe  given  in  the  following. 

For  ri2<yi2<  ti2,  that  is,  between  the  body  and  the  trailing  edge, 
it  is  seen  from  equations  (39b)  and  ( 3l)  that 


w 


-  V^a 


2yi 


■  2  “I 
1 


bx  -yx 
si2-yi2 


ACpi 


2V 


oa  ^s^-rj,2  (EjL-kj^1  K]_) 


>  (67a) 
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For  ti2<y12<si2,  that  is,  on  the  ving,  it  is  seen  from  equations  (39a) 
and  (29)  that 

1 


w  =  -VQa 


A9x  =  2V0a  A/s12-r12 


E(ki,tx)  -  kx  F(k1,^1) 


(67b) 


Finally,  for  yx2<  r i2,  that  is,  on  the  body,  equations  ( 38a)  and  (47) 
give 


w  =  -VQa 


Acpx  =  2V0a  y^*2 


t2 

kl  K! 


+  2V0a 


Ei-E(k1,t3)  +  kj_'  F(kx,t3) 

/(si2-y12)(r12-y12) 

J  . 


>  (67c) 


,  The  substitution  of  equations  (67a),  (67b),  and  (67c)  into 
equation  (66b)  yields  after  integration 


^2  “  -  MB12.r12)(Erk1,2K1)(E1,.k12K1’) 


(68) 


where  L/qa  is  given  in  equation  (65).  In  the  |  plane  equation  (68) 
can  be  written  in  the  dimensionless  form 


4qa 


2„  2 


8qas02 


“  (~T“)  (Eo-ko,2Ko)  (E0*  -ko2Ko* )  (69) 


which  for  rQ  =  0  agrees  with  the  results  of  reference  4.  Equation  (69) 
also  checks  with  the  result  obtained  for  the  drag  by  the  method,  pre¬ 
sented  in  the  preceding  section  on  section  drag,  based  on  the  calculation 
of  the  suction  force  along  the  wing  leading  edge. 
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(U) 


Sketch  (u)  shows  the 
total  drag  on  various  wing- 
body  combinations. 

Chord  loading.-  In 
order  to  find  the  center  of 
pressure  and  pitching  moment, 
it  is  convenient  to  find 
first  the  chord  loading, 
which  we  will  define  as  the 
value  of  /( Ap/qa)dy  where 
the  integration  is  carried 
over  the  wing  and  body.  The 
chord  loading  can  also  be 
obtained  by  evaluating  the 
expression  d(L/qa)/dx 
since  the  latter  term  is 
equal  to  /(Ap/qa)dy.  As 
a  check,  both  methods  were 
used  to  derive  the  following 
expressions. 

For  the  part  of  the 
chord  loading  contributed 
by  the  wing  it  can  be  shown 
that 


% 


A 


Mso4-r04) 


m  it-  arc  sin 


>  r0/m<x<c0 

02+r02/ 


(70a) 
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(70b) 


For  the  part  contributed  by  the  body  behind  the  leading-edge  fuselage 
juncture  it  can  be  shown  that 

p  a  d  (Lb^  4(s 04-r 04)  '  .  2s0r0  /  ^  ^ 

J  \<xU  ax  qa  So3  s02+ro2  - 


(71a) 


/Ap\  dy  =  _d_  _  4(s04-r04) 

\qai  dx  qa  s03 


V  qayb  "  dx  qa 
o 


m 


Eo  ^ 

K  p  -  E  p 

2  2 


(71b) 


The  total  chord  loading  can  be  obtained  by  combining  equations  (70) 
and  (71).  There  results  the  expressions 


for  r0/m<x<c0 


(72a) 


(72b) 
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Sketch  (v)  shows  the  variation  of  the  chord  loading  with  x 
for  r0/mc0  equal  to  0.316. 


O  4  .8  1.2  1.6 

(v)  x 


Center  of  pressure.-  The  results  of  the  last  section  can  he  used  to 
determine  the  center  of  pressure  xc>p>.  The  value  of  xc.p.  is  given 
hy  the  equation 


xc.p. 


L 


(73) 


which  excludes  the  loading  on  the  nose.  By  means  of  equation  (73) 
sketch  (w)  was  constructed. 
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II  -  ADDITION  OF  A  HORIZONTAL  TAIL 


It  is  possible  to  use  the  calculations  given  in  the  first  part  of 
this  report  to  find  the  forces  and  moments  induced  on  a  horizontal  tail 
by  the  presence  of  the  wing  and  body.  The  same  assumptions  that  were 
used  for  the  solution  of  the  load  distribution  over  the  wing  and  body 
will  be  made  here.  Hence,  the  results  will  be  principally  valid  for  air¬ 
planes  having  highly  swept  wings  and  tails  or  flying  at  Mach  numbers 
approaching  1. 

In  addition  to  the  basic  assumptions  by  which  slender  wing  theory 
is  defined,  however,  some  additional  assumptions  must  be  made  concerning 
the  behavior  of  the  vortex  sheet  trailing  behind  the  wing  and  passing  by 
the  tail.  Actually  these  trailing  vortices  provide  the  only  means  by 
which  the  wings  can  signal  their  presence  to  the  tail,  and  except  for 
them  the  slender  wing  theory  analysis  of  the  tail  effectiveness  would  be 
identical  to  that  described  in  part  I  for  the  wing.  Only  two  types  of 
trailing  vortex  patterns  will  be  investigated.  One  composed  of  a  flat 
vortex  sheet  situated  entirely  in  the  z  =  0  plane  (the  plane  of  the 
wing) ,  and  the  other  composed  of  two  completely  rolled  up  point  vortices 
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situated  symmetrically  with  respect  to  the  y  =  0  plane,  located  a 
distance  h  above  the  z  =  0  plane  and  a  distance  a  from  plane  of 
symmetry.  These  patterns  represent  the  two  extremes  of  the  actual  phys¬ 
ical  behavior  of  a  trailing  vortex  wake.  It  is  to  be  expected  that  the 
sheet  is  more  representative  of  the  true  wake  when  the  tail  is  located 
only  a  short  distance  behind  the  wing.  On  the  other  hand,  the  two  point 
vortices  should  be  valid  for  tails  located  a  large  distance  behind  the 
wing.  An  indication  of  the  magnitude  of  the  distances  at  which  the  two 
assumptions  are  accurate  can  be  obtained  in  reference  10. 


Method  of  Solution  and  Boundary  Conditions 


The  partial  differential  equation  that  governs  the  flow  in  the 
vicinity  of  the  tail  is,  of  course,  identical  to  the  one  studied  in  the 
first  part  of  this  report,  namely,  Laplace* s  equation  applied  to  a  yz 
plane  (equation  (2)).  In  fact,  the  general  discussion  of  boundary  con¬ 
ditions  and  forms  of  solution  given  in  part  I  still  applies  here.  Hence, 
the  Joukowski  transformation  can  again  be  used,  the  £  plane  having  the 
same  relation  to  the  plane  as  before  and  the  integral  relationship 

given  as  equation  (15)  still  applying. 

The  only  mathematical  difference  between  the  study  of  the  wing  and 
tail  can  be  seen  at  once  in  the  application  of  equation  (15).  In  the 
case  of  a  flat  wing,  the  vertical  induced  velocity  w  in  equation  (15) 
was  known  to  be  a  constant  over  the  region  occupied  by  the  wing  plan 
form.  In  the  case  of  a  flat  horizontal  tail,  on  the  other  hand,  the 
value  of  w  over  the  region  occupied  by  the  tail  plan  form  is  composed 
of  two  parts:  one,  the  constant  value  fixed  by  the  inclination  of  the 
surface  to  the  free  stream,  and  the  other,  a  distribution  that  is  just 
equal  and  opposite  to  the  vertical  velocity  induced  over  the  region  by 
the  vortices  trailing  from  the  wing.  Effectively,  therefore,  the  analy¬ 
sis  of  a  horizontal  tail  is  the  same  as  that  for  a  wing  with  a  given 
variation  of  twist  and  camber. 

The  additional  notation  necessary  for  the  description  of  the 
pertinent  tail  parameters  is  shown  in  sketch  (x) .  The  distance  from 
the  x  axis  to  the  tail  leading  edge  is  represented  by  a,  and  the 
slope  of  the  tail  leading  edge  is  designated  by  n.  In  this  report, 
only  triangular  tail  shapes  will  be  considered;  however,  more  complicated 
shapes  could  be  analyzed  by  the  method  presented. 
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Solution  for  Trailing  Vortex  Sheet 


Since  the  vortex  sheet  from  the  wing  is  assumed  to  lie  entirely  in 
the  z1  =  0  plane,  and  since  the  outer  extremities  of  this  sheet  are 
at  ±s^  (see  sketch  (y)),  the  study  of  this  case  can  commence  with  the 
inversion  of  the  integral  equation  (15)  of  part  I. 


wa(yi)  =  -  — 

2rt 


1  Av^ya) 


yx-y2 


For  t12<y12<s12,  the  value  of  Avx(y2)  is  given  by  equation  (25a)  and 
for  o12<y12<t12. 


Av1(y2)  =  0 


(75) 
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,z  Physical  plane 

Wing  vortex  wake 


,zt  Transformed  plane 


Wing  vortex  wake 


Substitute  these  two  values  of  Avx  into  equation  (7*0  and  apply  the 
boundary  condition  that  wj_  equals  -VQa  in  the  interval  0<y12<a12 
Then,  assuming  tj_2>aj.2,  that  is,  the  vortex  sheet  from  the  wing  does 
not  cross  the  tail10  (the  condition  shown  in  sketch  (y)),  there  is 
obtained,  after  inversion  (see  appendix  B)  and  some  manipulation,  the 
value  of  Av1  on  the  tail.  Thus 

for  0<y12<a12  and  for  tx2  >ax2 


Av1(y1)  =  - 


2V0ayi 


2V0ayi 


V-yi2  2 

zi 


dT)  /(n-'t12)(Ti-a12) 

y^-n  J  (s12-T])(T)-r12) 


This  solution  for  can  now  be  used  to  determine  the  aerodynamic 

characteristics  of  the  tail  in  the  presence  of  the  vortex  sheet  trailing 
back  from  the  wing. 


}This  assumption  applies  to  all  subsequent  analysis  of  the  tail  and 
vortex  sheet  combinations. 
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Span  loading . -  The  spanwise  variation  of  circulation  generated  by 
the  tail  surface  is  given  by  the  expression 


Acpt  = 


Av(y2)  dy. 


(77) 


wa 

where  Av  is  given  in  the  |x  plane  by  equation  (7 6).  On  the  portion 
of  the  plane  that  is  covered  by  the  tail  surface  (i.e., 

for  r12<y12<012)  this  yields 


2V0a 


n-*!3 


(s12-T])(n-r12) 


arc  tan 


'0i2-Ji3 


n-oi 


drj 


(78) 


To  determine  the  span  loading  on  the  body  it  is  necessary  to  subtract 
the  value  of  at  the  tail-leading-edge  fuselage  juncture.  This 

value  is  obtained  from  equation  (78)  by  placing  =  r±.  The  span  load¬ 
ing  on  the  body  (i.e.,  for  0<y12<r12)  is  then  given  by  the  formula 


2V0a 


n-ti2 

-Tl)(Tj-r12) 


arc  tan 


■  /a12-y1£ 

/ - —  dt) 

J  T)-Ol 


(79) 


Equations  (78)  and  (79)  have  been  transformed  to  the  |  plane  and  the 
results  are  shown  in  figure  2(a).  In  this  and  in  all  following  numeri¬ 
cal  examples  the  wing  will  be  fixed  as  the  special  type  studied  in 
part  I  having  the  measures  s0/mc0  =  1.7,  t0/mcQ  =  1.091  and 
r0/mc0  =  0.316. 
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Chordwise  load  distribution  on  the  tail.-  The  distribution  of  load 
over  the  tail  surface  can  be  calculated  from  the  equation 


/ Ap  \  2  clAcpt  _  2  cWPit  doi 

\  qa^  VQa  8x  VQa  8a  x  dx 


(80) 


where  the  value  of  lt  determined  from  equation  ( 78) ,  is  given 

by  the  expression 


8A9i+  2VoOUJ,  2VQ0WJ, 

—  —  1  ■  +  „  —  

S(Ji  *  Jo? -y? 


(il-ti2)  drj 


V 


y(i Sl2-n)  (  11-r!2)  (  TJ-O^)  (  TJ-tx2) 


(81) 


By  means  of  the  substitutions 


(g1g-rig)(si2-t1g) 

(s1g-a1g)(t1g-r1g) 


sngu 


(si2-gig)(n-ti2) 

(si2-t12)(T]-012) 


(82) 


equation  (8l)  can  be  reduced  to 

8Aq>lt  = 

^1 

where 

+9 

Hence,  the  loading  on  the  tail  < 


2V0cw1A0(k5,t9) 


be  written  in  the  closed  form 


(83) 


(84) 


and  since 


ik^ApUs,1^)  do1 

V° i2_yi2  ** 


(85) 


« 


* 
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^7 


the  final  expression  obtained  by  transforming  equation  (85)  to  the 
|  plane  is 

(%)  -  h2l|^A0(k4,»a)  (86) 

\qa/t  aJ 

where  k4  and  ife  are  the  transformed  values  of  k5  and  +9  and  are 
given  in  the  table  of  symbols. 

Total  lift  on  the  tail.-  The  total  lift  on  the  tail  can  be  evaluated 
by  use  of  the  equation 


Lt 

qa 


a 


o 


A<P(y)  dy 


which  becomes  in  the  £1  plane. 


Lt 

2qa 


1 

VQa 


ox 


y1Acp1(yi)dy1 


/ 


yi2-ri2 


(87) 


(88) 


The  value  of  ATi  is  given  by  the  equations  (78)  and  (79)*  Substitute 
these  expressions  into  equation  (88)  and,  after  integration,  there 
results 


=  it(t12-r12-s12+a12)  AQ(k5,+9)  + 

2qa 


2y(S^i=)(t^) 


t1=-ol2  , 

E=  -  (  — - —  1 


s  x2-Ox2/  J 


(89) 


In  the  £  plane  this  can  be  written 


Lt. 

2qa 


=  —  it 


( s02-t02) ( s02t02-r04)  (a02-r02) 


so2to2 


*o2 


A-o(^4>'lf8^  + 


1 2 

- - - - -  J  (so2“0o2)(  s02002“r04)  “  4f2 

t0S0  V  \  k0*  / 


(90) 


A  plot  of  this  equation  is  shown  in  figure  3(a)  for  r0/mc0  =  O.316. 
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Drag.-  The  drag  of  the  tail  can  he  calculated  from  the  equation 


Dt  =  aLt  +  Ft  (91) 

where  L-t  is  given  by  equation  (90)  and  Ft  is  the  suction  force  at 
the  leading  edge.  In  a  manner  similar  to  that  given  in  the  first  part 
of  this  report  under  the  subheading  "Section  drag,"  this  force  is 
obtained  from  the  equation 


(92) 


As  before 


G  =  lim  B  un(xn>yn) 
xn-S>  o  n  V0 


(93) 


the  values  of  un  and  ^  being,  respectively,  the  normal  velocity  to 
and  the  normal  distance  from  the  leading  edge.  Substituting  these  values 
into  equation  (93)  gives 


G2  = 


a 


2vC7  v  y3 


y  -*< 


An2(k 


ts) 


4**8 


(94) 


and  the  expression  for  the  suction  force  can  he  written 


Ft  =  -  V02a2 pit 


Aos(k4^8)3y  (95) 


Hence,  the  total  drag  becomes 


-  (96) 

ct  q  aq  J  \  a3  / 

ro 

This  equation  is  plotted  in  figure  4(a)  for  r0/mc0  =  0.316. 


Chord  loading.-  A  closed  formula  for  the  loading  can  be  obtained  by 
carrying  out  the  integration  /  (Ap/qa)tdy  over  the  tail  or  by  evaluating 
t^(4t/<la')/dx.  However,  since  the  term  Ao(k4, 1lra)  in  the  expression 
for  Ap/qa  does  not  involve  y,  it  is  easier  to  evaluate  /(Ap/qa)tdy. 
Thus,  using  equation  (85),  the  chord  loading  in  the  plane  is  given  by 
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(97a) 


In  the  i  plane  this  becomes 

y-w(^)Ao(^*J  '  (97b) 

o 

(This  checks  with  the  result  found  by  differentiating  Lt/qa  with 
respect  to  x) . 


Center  of  pressure.-  The  center  of  pressure  can  be  calculated  by 
means  of  the  formula 


xc  .p . 


Lt 


(98) 


where  dLt/dx  is  the  result  just  obtained  in  the  previous  section 
and  Lt  is  given  by  equation  (90).  Placing  these  values  into 
equation  (98)  yields 


•P-  ■  TT  /  Ao<^,*a)  to 


(99) 


A  graph  of  this  result  is  shown  in  figure  5(a)  for  r0/mcQ  =  0.316. 


Solution  for  Rolled  Up  Vortices 


As  was  pointed  out  in  the  preceding  section  where  the  method  of 
solution  was  discussed,  the  manner  in  which  the  present  problem  will  be 
attacked  is  as  follows:  First,  the  velocities  induced  at  the  surface  of 
the  wing  and  body  by  two  point  vortices  located  somewhere  in  space  will 
be  calculated;  second,  a  solution  will  be  formulated  (by  methods  identi¬ 
cal  to  those  used  in  part  I  of  this  report)  that  will  just  cancel  the 
vortex  induced  velocity  component  normal  to  the  surface  of  the  wing  or 
body;  third,  an  additional  solution  will  be  formulated  that  will  fit  the 
boundary  conditions  prescribed  for  the  tail  surface  (in  this  report  only 
a  flat-plate  tail  surface  will  be  considered). 
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As  usual,  it  is  simpler  to  work  with  the  transformed  plane 

than  with  the  physical  £  plane.  Hence,  again  the  Joukowski  transfor¬ 
mation  will  be  applied  to  the  field  equation  and  boundary  conditions  at 
the  outset  of  the  problem.  See  sketch  (z). 


The  velocity  potential 
at  a  point  (yi,zi)  in  the 
| 1  plane  induced  by  a  pair  of 
point  vortices  located  at 
yi  =  ±ai  and  z*  =  hx  is 
given  by  the  equation 

cp  =  £w  r&1  (z1-h1)dyg 

2 *  J  (yi-y2)2+(zi-hi)2 

"al 

(100) 


Transformed  plane 


Vortex  location— r 

1/ 


where  is  the  strength  of 

the  circulation  carried  by  the 
wing  panels  and  trailing  back 
from  the  wing  tips.  The  value 
of  rw/v0<xmc0  that  corresponds 
to  the  swept-back  wings  studied 
in  part  I  of  this  report  is 
given  by  equation  (32a),  thus 


^w 

V0amc0 


=  2 


The  values  of  ax  and  h^  depend,  of  course,  on  the  span  of  the  wing. 

In  order  to  compare  the  following  results  for  the  rolled  up  vortices  with 
those  obtained  for  the  sheet  vortices  in  the  preceding  section,  the 
numerical  results  presented  in  the  succeeding  examples  will  be  for 
s0/mc0  =  1.7,  t0/mc0  =  1.091,  and  r0/mc0  =  0.316.  This  particular 
choice  of  parameters  fixes  the  span-loading  curve  for  the  wing  to  be  that 
shown  in  sketch  (o).  A  reasonable  choice  for  the  value  of  a  can  be 
calculated  by  replacing  the  figure  in  sketch  (o)  bounded  by  the  lines 
y  =  r0 ,  T  =0  and  the  curve  for  T  /v0amc0  by  a  rectangle  with  the  same 
height  and  area.  The  value  of  a  is  then  given  by  the  sum  of  the  base 
length  of  this  rectangle  and  the  quantity  r0/mc0.  This  procedure  was 
carried  out  and  the  result  a/mc0  =  1.5^5  was  obtained.  In  order  to 
obtain  a  more  complete  picture  of  the  effect  of  the  point  vortices  on 
the  aerodynamic  characteristics  of  the  wing,  four  different  locations 
were  chosen  for  the  positions  of  the  vortices,  two  in  the  z  =  0  plane 
for  values  of  a/mc0  equal  to  1-5^5  and-  1.3>  and  two  at  a  height 
h/mc0  =0.3  above  the  z  =  0  plane  for  the  same  two  values  of  a/mc0. 
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From  equation  (lOO)  it  can  easily  be  shown  that  the  value  of  the 
vertical  induced  velocity  in  the  z±  =  0  plane  is 


\  _  iv 

yi-ax 

yi+ai 

^zi4,=0  2* 

Lhi^Cyi-ai)2 

ha2+(yi+a1)2  J 

(101) 


The  particular  inversion  of  equation  (15)  that  fits  the  present  boundary 
conditions  ( see  appendix  B)  can  be  written 


A Viiyi)  =  — ------  p01  J  ai2-y22  dy2  (102) 

±ax  yi‘ys 

Place  the  value  of  wv  given  by  equation  (lOl)  into  equation  (102)  and 
add  the  condition  that  the  tail  be  a  flat -plate  lifting  surface  at  an 
angle  of  attack  a,  and  there  results 


Av-^y-L) 


J  hi2 -Ci2  yb^-ai2 

—————— 

yi-bj.  yi-bi 

where 

bi  =  ai  +  ihx 

bi  =  ai  -  ihi  (104) 

and  where  the  radicals  are  defined  uniquely  if  the  complex  plane  is  cut 

along  the  real  axis  between  -  °°  and  o1,  (For  example,  can 

be  set  equal  to  p1e^cPi  where  cpx  must  lie  between  -it  and  it.) 
Although  the  above  expression  for  Av  can  be  put  in  real  terms  by 
applying  the  transformations 

a  1  =  <?!  cosh  y1  cos 

| hi  |  =  ox  sinh  7±  sin  ,  (105) 

it  is  easier  in  deriving  subsequent  quantities  from  Avx  to  use 
equation  (103)  first  and  to  make  the  transformation  afterwards. 


-2Yqclj1 


r  w 


J  a12,*yi2  <*1 


2-yi2 


/bis-Oi« 


— —  + 


yi+bi 


yi+bi 
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The  terms  sinh  yl9  sinh  sin  co^  and  sin  uj2  are  all  defined  in  the 
table  of  symbols.  Equations  (107a)  and  (l07b)  can,  of  course,  be  trans¬ 
formed  to  the  physical  plane  by  transforming  each  symbol  therein  from 
the  lx  system  to  the  I  system. 

The  variation  of  the  tail  span  loading,  as  given  by  equations  (107a) 
and  (107b),  is  shown  in  figures  2(a)  and  2(b)  for  the  various  vortex 
positions  discussed. 

Chordwise  load  distribution.-  The  loading  on  the  tail  can  be 
calculated  from  the  relation 


/Ap\  2  2  ^Pit  dai 

Vqa/j.  V0a  dx  VGa  dx 

The  value  of  dATit/cto^  follows  from  equations  ( 107)  j  thus 


(108) 


d&Pit 


2V0cus1 


Jvi2 -Jiz 


1  + 


Pw 


Hence, 


(-) 

.  _ 

rw  sinh  7a  cos 

dai 

\  qaX 

t  Ta^-yi2 

V0a«a1(sinh27i+sin2u)1)  _ 

dx 

(109) 


(no) 


Equation  (110)  when  transformed  to  the  |  plane  gives  the  surface  load 
distribution  over  the  tail  due  to  the  presence  of  the  two  point  vortices 
as  well  as  the  inclination  of  the  tail  to  the  free  stream. 


Total  lift.-  The  total  lift  can  be  obtained  by  integrating  the  span 
loading.  Thus,  if  Lt  represents  the  lift  on  the  tail, 


r 

qa  VQa  J 


Carrying  out  this  integration  yields 


(HI) 


h  *£,  2P 

-  =  j^o-,2-^2)  -  — -  (r,  sinh  7~  cos  to„-an 

2qa  1  1  VQa  1  2  2  1 


sinh  /j^cos  w^) 


(112) 
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This  result  was  transformed  and  plotted  in  figures  3(a)  and  3(h)  for  the 
values  discussed. 


Total  drag.-  The  total  drag  is  given  by  the  relation 


Dt  =  aLt  +  Ft 


(113) 


where  Ft  is  the  suction  force  on  the  tail  leading  edge.  As  in  previous 
steps  (see,  e.g.,  equation  (92))  the  calculation  of  Ft  depends  on  the 
evaluation  of  the  function  G.  In  this  case  G2  is 


h-  — 

L  2V0ait 


2  (114) 


The  final  result,  in  the 


£1 


plane ,  for  the  total 


drag  can  he  written 


Dt_ 

ga2 


it(ai2-ri2)  +  i- 
2n 


( a x2-h xa-a  1«)s+4a12h x2  +( a x2  -h x2  -a x2 )  J (a12-h12-cT12)2+4a12h12 
( a12-h12  _ri2  )2+4a12h12+(  ax2  -hj2  -rx2 )  J~(  a12-h12-r12)2+4a12h1a 


(115) 


A  plot  of  the  drag  is  given  in  figures  4(a)  and  4(h). 

Chord  loading.-  As  before,  in  the  development  of  equations  (97a) 
and  (97b),  for  example,  the  chord  loading  can  he  calculated  either  by 
performing  the  integration  /(Ap/qa)-^  dy,  or  by  differentiating  with 
respect  to  x  the  total  lift  (Lt/qa) .  These  two  different  approaches 
serve  to  check  each  other  and  both  lead  to  the  same  result,  namely, 


2 


dy  =  4^! 


dcfi 

dx 


1  - 


rv  sinh  7j  00s 


V0aa1it(  sinh271+sin2aj1) 


(116) 


Center  of  pressure.-  Results  for  the  center  of  pressure  xc.p. 
(where  xc.p.  =  -Mt/Lt)  are  shown  in  figures  5(a)  and  5(b) . 
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CONCLUDING  REMARKS 


When  an  airplane  is  slender  enough11  (in  the  longitudinal  sense)  or 
is  flying  close  enough  to  the  speed  of  sound,  the  mathematical  descrip¬ 
tion  of  its  attendant  flow  field  is  greatly  simplified  -  so  much  so,  in 
fact,,  that  the  analysis  of  whole  wing -body-tail  combinations  is  feasible. 
This  simplification  comes  from  the  fact  that  the  induced  velocities  in 
each  lateral  plane  are  completely  independent  of  the  nature  of  the  air¬ 
plane  or  flow  field  behind  the  reference  plane  and  are  affected  by 
disturbances  ahead  only  through  the  presence  of  free  vortices  trailing 
downstream  from  the  lifting  elements.  In  the  case  of  a  tail,  these  free 
vortices  stream  back  from  the  wing  trailing  edge. 

In  this  report,  special  wing  plan  forms  were  studied:  special  in 
that  they  produced  flat  span-loading  curves  between  the  wing  tips  and 
fuselage.  For  such  wings,  the  free  trailing  vortices  were  concentrated 
entirely  in  the  region  directly  behind  the  wing  tips.  In  general,  the 
trailing  vorticity  would  be  concentrated  predominately  in  this  region. 

The  behavior  of  this  trailing  vortex  system  is  bounded  by  the  behavior 
of  two  extreme  models:  a  vortex  sheet  lying  everywhere  in  the  plane  of 
the  tail,  and  two  laterally  symmetric  point  vortices  lying  in  or  above 
the  plane  of  the  tail.  Each  of  these  models  was  examined. 

One  point  vortex  was  placed  in  the  plane  of  the  tail  at  a  distance 
from  the  fuselage  in  the  spanwise  direction  determined  by  replacing  the 
wing-span  loading  curve  by  a  rectangle  of  the  same  height.  As  shown  by 
figures  2  through  5*  the  results  for  this  point  vortex  were  not  signifi¬ 
cantly  different  from  those  for  the  vortex  sheet.  In  either  case,  the 
presence  of  the  trailing  wing  vortices  reduced  by  about  40  percent  the 
effectiveness  of  the  triangular  tail  surface  in  producing  lift  for  the 
range  of  tail  spans  and  body  diameters  considered.  For  the  same  condi¬ 
tions  the  tail  drag  was  reduced  only  18  percent. 

For  the  particular  locations  chosen  for  the  point  vortices,  it  was 
found  that  both  the  lift  and  drag  decreased  as  the  vortices  moved  closer 
to  the  tail.  On  a  percentage  basis  the  decrease  was  roughly  the  same. 


11The  assumptions  underlying  slender  wing  theory  are  obviously  violated 
along  lines  such  as  the  leading  edge,  x  =  mc0>  and  the  Mach  wave  from 
the  trailing-edge -fuselage  juncture,  x  =  c0.  Along  these  lines  the 
pressure  gradient  is  discontinuous  and  (Mq2-1)  cpxx  is  not  bounded. 
Similar  situations  appear  repeatedly  in  the  linearized  analysis  of 
aerodynamic  flow  phenomena  and  in  each  case  agreement  with  experi¬ 
mental  results  cannot  be  anticipated. 
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The  position  of  the  center  of  pressure  on  the  triangular  tail  was 
insensitive  to  the  presence  of  the  wing  vortex  system  regardless  of  the 
vortex  pattern  chosen.  In  the  extreme  case,  when  the  point  vortices  were 
nearest  the  tail,  the  location  of  the  tail  center  of  pressure  with  refer¬ 
ence  to  the  tail  apex  as  5  percent  forward  of  the  position  obtained  when 
the  wing  was  absent. 


Ames  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Moffett  Field,  Calif.,  Aug.  20,  1951 
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APPENDIX  A 

LIST  OF  IMPORTANT  SYMBOLS 


c0 

cd 

cz 

d 

D 


horizontal  distance  from  y=0  plane  to  vortex 
(See  sketch  (z) . ) 

ax  +  ihi 

-  ihjL 

local  chord 

characteristic  chord 
(See  sketch  (a),) 

section  drag  coefficient  [  — 

\  qc 

section  lift  coefficient  (  — 

\  qc 

section  drag  force 
drag  force 

complete  elliptic  integral  of  the  second  kind 


E(k,+)  incomplete  elliptic  integral  of  the  second  kind 


F  suction  force  at  leading  edge  of  lifting  surface 

F(k,+)  incomplete  elliptic  integral  of  the  first  kind 
*  dt 


(IJW-) 

( rti-k2t2~ 

LndU/l^- 


V(l-k2t2)(l-t2) 


h 

Kb 


k 

k» 


vertical  distance  from  z=0  plane  to  vortex 
(See  sketch  (z).) 

complete  elliptic  integral  of  the  first  kind 

(  r1 .  dt 

\Jo  /( l-k02)(l-t2) 
modulus  of  elliptic  integrals 

y  i-k2 


58 


NACA  TN  255k 


ko 

kx 

k2 

k3 

k4 

k5 

Z 


si2-tx2 


</  (sgt2-r04)(s2-t2) 
t(s2-r02) 

ko  in  the  lx  plane  ,  _ _ 

\  V  si2-ri2 

2roy  (s2tg-r04)(s2-tg) 

(s2-r02)(t2+r02) 

/k  rx*/si2-ti2  ^ 
k2  in  the  |x  plane  (  ^=j  J 

(g2-^2)^2-!2)^^*) 

(t2-r02)<y  (s2-a2)(s2CT2-r04) 

(<Jig-ni2)(si2-t12)~ 

T) 


k4  in  the 
section  lift  force 


*  .  /  y (^i2-ri2)(s12-t12)  "\ 

5i  piane  ; 


L 

m 


M 

Mo 

P 

4P 

ol 


s 


s0 


S 

t 

*0 


lift  force 

slope  of  wing  leading  edge 
(See  sketch  (a),) 

pitching  moment,  positive  when  tail  is  forced  down 
free-stream  Mach  number 
local  static  pressure 


loading  coefficient 


PZ  ~Pu  ^ 

J 


free-stream  dynamic  pressure 


radius  of  body 

(See  sketch  (a).) 

distance  from  x  axis  to  wing  leading  edge 
(See  sketch  (a).) 

maximum  value  of  s 
(See  sketch  (a).) 

wing  area 

distance  from  x  axis  to  wing  trailing  edge 
(See  sketch  (a).) 


* 


M 


l 


maximum  value  of  t 
(See  sketch  (a).) 
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u,v,v  perturbation  velocity  components  in  the  x,y,z  directions, 
respectively 


Au,Av,Aw  jimp  in  velocities  across  z=0  plane,  u,.-u7,  v..-v7,  , 

respectively 

vr  radial  component  of  perturbation  velocity  in  a  y,z  plane 

V0  free-stream  velocity 


x,y,z 

xc.p. 


a 

P 

sinh  7 
sinh  r2 

r 

n 

9 

AoM) 

H* 


I 

P 

Po 

a 


A<P 

(Acp)T>E_ 


Cartesian  coordinates 

distance  to  center  of  pressure 
(See  sketch  (w)  and  fig.  5.) 

angle  of  attack  of  airplane 

•J  Mo2"! 

Ihxl 


crjLsina)! 


rx  sin 


total  circulation  about  wing  section 


polar  angle  in  y,z  plane 

|  [KE(k',Y)  +  (E-K)P(k«,t)] 

slope  of  tail  leading  edge 
(See  sketch  (x).) 

complex  variable  (y+iz) 

polar  distance  in  y,  z  plane  ( */  y2+z^) 

free-stream  density 

distance  from  x  axis  to  leading  edge  of  tail 
(See  sketch  (x).) 

perturbation  velocity  potential 

jump  in  <p  across  z=0  plane  (cPu-q^) 

value  of  A 9  at  trailing  edge 
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t 

*0 

+. 2 

*3 


argument  of  elliptic  integrals 

«/(s^a-rQ4)(sg-y2)t^ 

^  (  s^t  2_r0  4)  (  s  2-t  y  2 

+0  in  gx  plane  f  ) 

V  </  Si2-!/  / 

2t</rQg-yg _ 

a/  (t2+r02)2  -  l^t2 

♦2  in  ii  Plane  ( 

\  Jti 2-yis 


+4 

*6 

+7 

*8 
*9 
.  2 

sin  a)! 
sin2  u)2 


s(t  +rQ  ) 
t(s2+r02) 

t4  in 


2  2 
s  -r0 

2  2 

s  +rQ 


in  lx 


s(-t^-rQg) 

t(^-r02) 

ts  in  li  plane 


(  -j  tia-ng  \ 

n  V  si2  -if2  / 


[  0ig-a12-h1g+«/(a1g+h1g-g1g)2  +  ko^h^) 


2ai2 

[rig-ai2-hi2+</(aig+hig-r1g)g  +  it-rx2^2] 

2rig 


V 


Subscripts 

1  complex  plane  resulting  from  the  application  of  the  Joukowski 

transformation  to  the  physical  plane 

12  boundary  conditions  given  by  equations  (12) 

t3  boundary  conditions  given  by  equations  (13) 
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a 

b 

Z 

n 

t 

u 

v 


wake 

body- 

lower  surface  of  z=0  plane 
component  normal  to  leading  edge 
tail 

upper  surface  of  z=0  plane 
rolled  up  yortex 


w 


wing 
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APPENDIX  B 

INVERSION  OF  AN  INTEGRAL  EQUATION 


The  integral  equation 


f(ni) 


1  r  Av(T]a)dng 

2it  Ja  TJ  i_TJ  ^ 


(Bl) 


can  be  inverted  by  applying  operational  techniques, 
operator 


h(m)ani 

X-R  1 


Consider  the 


(B2) 


where  h(r|1)  is  a  function  to  be  chosen  later.  Operating  on  both  sides 
of  equation  (Bl)  yields 


Mu  3)  an 

X.-1U 


1 


AvhsMui) 

(r  1-TJ2)(X-T1  1) 


(B3) 


where  the  order  of  the  integral  and  differential  signs  on  the  right-hand 
side  of  the  equation  indicates  that  the  t\z  integration  is  to  be  per¬ 
formed  first. 

The  next  step  is  to  reverse  the  order  of  integration  in  the  double 
integral  term  in  equation  (B3).  Since  an  inherent  singularity  exists  in 
the  area  of  integration  at  the  point  T]i  =  tj2  =  A.,  however,  some  care 
must  be  used  in  order  to  obtain  this  reversal.  Designating  by  R(\)  the 
difference  between  the  term  taken  first  with  one  order  of  integration  and 
then  with  another,  thus 

K(x)  -  r\jb Avf^hM  _r\z  ATMhOu) 

Ja  Ja  (m-RaUX-Ti  1)  Ja  Ja  (n  i-na)(x-iu) 

(B4) 

R(x)  can  be  evaluated  by  isolating  the  singularity  and  studying  the 
difference  only  in  its  vicinity.  Hence, 
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R(x)  =  lim 
£->0 


X+e  X+e 

P  drji  P  dt]2 


Av(Tig)h(ru) 

(Tll-'rl2)(x-T)l) 


X+e  X+e 
f  dT]2  P  din 


Av(ng)h(T|l) 

(t!i-t12)(x-tu) 


Using  the  transformations  ^  =  X-€xx  and  r|2  =  X-ex2,  equation  (B5) 
can  be  reduced  to  the  form 


R(X)  =  Av(x)h(x)  P  dxiT  -  P  dx2P  - ^ — 

Lli  J-1  Xl(x2-Xi)  Jj_  Z-i  x1(x2-x1) 


=  2Av(x)h (x)/°  ^  Zn 

J0  x  1  1+Xl 

=  -rt2Av(x)h(\)  1 

Substitute  equations  (B6)  and  (BU)  into  equation  (B3)  and  there 
results 


b  h(Tu)f(Tu)dni  _  _  1  r  _„2a rb  h(ni)Mi 


f  Av(x)h(x)  +  /  av^W  -iA-'ii  ^11. 

2ir  L  ,  Ja  *4  ( Tii  -ri2 )  (x-tii  ) 


Av(x) 


2  r  1 


:h(x)  L  at 


-AvtyXlna/*  il(’ll)a11 - +  /*h(n)f lnAau  l  (B7) 

<~4  J&  ( 'Hi  -Tl2  H  X-T)i  )  4,  X-T)l  J 


If  h(ni)  =  fib  -’rl1)(T]1-a)  equation  (B7)  becomes 


Av(x)  =  - -==.:2,,,,:,  -  P  Av(  T|2  )  dTJ2  +  f  ^(b  "ll  )  ~a')  dTJ  x 

*a/( b-x)(x-a)  -  2Ja  Ja  X-m 


Since,  hovever,  equation  (b8)  contains  both  Av  and  its  integral,  it 
does  not  represent  a  unique  solution  for  Av.  In  order  to  obtain  a 
unique  solution,  some  additional  condition  must  be  imposed.  If  this 
condition  is  that  Av(a)  =  0,  then  the  proper  choice  of  h(r|1)  in  equa¬ 
tion  (B7)  is 


h(n1)  =  /IH 
V  b-Tl! 

which  leads  immediately  to  the  unique  result 


A  v(x) 
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f°  f(Ti  i3dn  i 
h  -X  Ja  X-T)  x 
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TABLE  I.-  POSITION  OF  TRAILING  EDGE 


r0/mc0  = 

0 

0.100 

0.224 

0.316 

0.447 

0.500 

0.600 

0.700 

0.800 

0.900 

x/  c0 

y/mc0 

y/mc0 

y/mc0 

y/mc0 

y/mco 

y/mc0 

y/mc0 

y/mco 

y/mc0 

y/mc0 

1.00 

0 

0.100 

0.224 

0.316 

0.447 

0.500 

0.600 

0.700 

0.800 

0.900 

1.05 

.168 

.213 

•  323 

.413 

•  542 

•  593 

.690 

.791 

.883 

•  970 

1.10 

.256 

.290 

.387 

.473 

.598 

.649 

.746 

.844 

.933 

1.020 

1.15 

.332 

.361 

.447 

.526 

.648 

.699 

.796 

.890 

.982 

1.067 

1.20 

.400 

.425 

.504 

.579 

.698 

.748 

.845 

•  933 

1.031 

1.113 

1.30 

.526 

.545 

.614 

.682 

.797 

.845 

•  941 

1.028 

1.126 

1.213 

1.40 

.644 

.661 

.721 

.786 

.896 

.941 

1.033 

1.123 

1.221 

1.312 

1.50 

.756 

.773 

.826 

.888 

.993 

1.036 

1.128 

1.217 

1.315 

1.410 

1.60 

.867 

'  .884 

.931 

•  990 

1.090 

1.130 

1.222 

1.312 

1.412 

1.506 

1.80 

1.083 

1.096 

1.140 

1.192 

1.283 

1.323 

1.413 

1.504 

1.606 

1.703 

2.00 

1.297 

1.308 

1.394 

1.393 

1.482 

1.522 

1.610 

1.700 

1.800 

1.800 

2.20 

1.507 

1.516 

1.555 

1.596 

1.676 

1.723 

1.803 

1.893 

1.996 

2.100 

2.60 

1.91 

1.92 

1.96 

2.00 

2.07 

2.12 

2.20 

2.29 

2.39 

2.50 

3-00 

2.32 

2.33 

2.37 

2.40 

2.47 

2.52 

2.60 

2.69 

2.79 

2.90 

3.50 

2.83 

2.84 

2.87 

2.91 

2.98 

3.02 

3.10 

3.19 

3.29 

3-39 

4.00 

3.33 

3-34 

3.38 

3.41 

3.49 

3.52 

3.60 

3.69 

3.78 

3.89 

s 


[ilw 


gjag 


— 

0.100 

-4- 

CVJ 

OJ 

# 

0 

VO 

H 

OO 

• 

O 

0.447 

0.500 

O 

o\ 

O 

O 

0.700 

0.800 

0.900 

y/mc0 

C/C0 

c/c0 

c/cq 

c/c0 

C/C0 

c/c0 

c/c0 

c/c0 

c/c0 

0 

1.000 

-  -  - 

_  .  _ 

_  _  _ 

_  _  _ 

mm  mm 

... 

mm 

.100 

.926 

0.900 

-  -  - 

-  -  - 

—  —  _ 

_  _  _ 

mm  mm  m. 

... 

... 

mm  mm  ^ 

•  22J 

.910 

.831 

0.776 

-  -  . 

-  -  - 

... 

mm  mm  . 

... 

mm  ~  mm 

mm  mm  mm 

.316 

.826 

.802 

.732 

0.684 

-  -  - 

—  -  — 

... 

... 

mm  mm  — 

mm  mm  mm 

.447 

•  790 

.772 

.703 

.633 

0.553 

-  .  _ 

... 

... 

... 

mm  mm  — 

.500 

•  777 

.761 

.695 

.627 

.522 

0.500 

... 

... 

... 

.600 

.760 

.746 

.685 

.620 

.504 

.462 

o.4oo 

... 

... 

... 

.700 

•748 

•  735 

.679 

.617 

.498 

.451 

.362 

0.300 

... 

... 

.8 

•  738 

.724 

.673 

.614 

.500 

.452 

•  355 

.262 

0.200 

... 

•  9 

.728 

•715 

.668 

.612 

•  503 

.459 

•  358 

.260 

.170 

0.100 

1.0 

.722 

•  709 

.665 

.610 

.505 

.464 

.365 

.269 

.170 

.082 

1.1 

.719 

.704 

.661 

.608 

.508 

.468 

.372 

.278 

.175 

.079 

1.2 

.708 

.697 

.658 

.607 

.512 

.472 

•  376 

.285 

.180 

.081 

1.3 

.703 

.692 

.653 

.606 

.515 

.475 

.381 

.290 

.185 

.087 

1.4 

.699 

.688 

.651 

.605 

.518 

.477 

.385 

.294 

.189 

.094 

1.6 

.692 

.682 

.645 

.603 

.522 

.480 

.392 

.299 

.196 

.098 

1.8 

.687 

.677 

.642 

.601 

.524 

.481 

•  396 

•  302 

.200 

.100 

2.0 

.682 

.672 

.639 

.599 

.525 

.482 

.398 

•  305 

.204 

.100 

2.5 

.674 

.663 

.631 

.595 

.522 

.484 

.403 

.311 

.212 

.101 

3.0 

.669 

.657 

.625 

.593 

.519 

.485 

.406 

•  315 

.218 

.102 

4.0 

.660 

.648 

.613 

.592 

.514 

. 

.485 

.410 

.318 

.225 

.102 
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Figure  /.  -  Relations  between  the  wing  area,  wing  span,  and 

body  diameter. 
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O  .2  4  .6  .8  1.0 

/s„ 


(b)  Constant  t0/s0 
Figure  I.  -  Concluded. 


'  Vcfi.mco  represents 
loading  developed 
behind  section  AA 


Variation  of  the  span  loading  on  the  tail 


(b)  Point  vortices. 
Figure  2.  -  Concluded 


Of* 


NACA  TN  255^ 


<r0/mc0 


2  ,  4  .6  .8  1.0 

eg  Fmc0 


(b)  Point  vortices. 

Figure  3.  -  Variation  of  the  lift  on  the  tail. 
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(a)  Vortex  sheet. 


(b)  Point  vortices. 

Figure  4.-  Variation  of  the  drag  on  the  fait. 
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